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We simulate 3-flavour lattice QCD with a small chemical potential fii for isospin, at temperatures close to the 
finite temperature transition. Using quark masses just above the critical mass for zero chemical potential, we 
determine the position of the transition from hadronic matter to a quark-gluon plasma as a function of We 
see evidence for a critical endpoint where the transition changes from a crossover to a first-order transition as 
fii is increased. We argue that QCD at finite and QCD at finite quark-number chemical potential \i should 
behave similarly in this region. 



1. Introduction 

QCD at finite temperature and quark-number 
chemical potential, (i, is of relevance to the 
physics of relativistic heavy-ion collisions. Dif- 
ferent parts of the phase diagram are accessible 
to the various heavy- ion accelerators — RHIC, 
CERN heavy-ion, AGS and SIS. 

For the region of small (i and temperatures 
close to the finite temperature transition, various 
techniques have been developed to circumvent the 
problems associated with the complex fermion de- 
terminant at finite fj, |1I2I3I4| . We have chosen 
to study the finite temperature transition with 
a finite chemical potential /i/ for isospin, which 
is equivalent to simulating QCD at finite fi with 
2/i = fij using only the magnitude of the determi- 
nant and ignoring the phase. Because the fermion 
phase is relatively well behaved at small /i, it can 
be argued that the positions of the transitions 
and probably their nature with and without this 
phase, should be the same. 

For 2 flavours we found no sign of a critical end- 
point for /ij < TOtt, the maximal range over which 
finite /ii and finite /i transitions could be expected 
to be related [5] . Here we describe our simulations 
for 3 flavours. Not only is this more physical, but 



*Talk presented by D. K. Sinclair, Lattice2004, Fermilab. 
t Supported in part by NSF grant NSF PHY03-04252. 
*Supported by DOE contract W-31-109-ENG-38. 



for m < m c (0), the finite temperature transition 
is first order at /// = 0. At m = to c (0) it becomes 
second-order, in the Ising universality class. For 
m > m c (0) the /i/ = transition is a crossover. It 
is expected that the critical endpoint moves con- 
tinuously from \ii = for to = to c (0) to finite ///, 
as to is increased. Thus we should be able to keep 
this endpoint as close to /ij = as we please, by 
choosing m to lie just above m c (0). The prelim- 
inary results of our simulations at such a mass, 
reported here, do show such a critical endpoint 
at a modest value of /ii. 

In section 2 we describe our simulations and 
report preliminary results. Our conclusions are 
given in section 3. 

2. 3-flavour QCD at finite fij and T 

The staggered fermion part of our lattice action 
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For these simulations, which are at /i/ < m w , we 
set the symmetry breaking parameter A = 0. 

We are performing simulations on 8 3 x 4 and 
12 3 x 4 lattices at m = 0.03, 0.035 and 0.04. 
[For N t = 4, m c (0) = 0.0331(12) ©.] We use 
standard hybrid molecular-dynamics simulations. 
Since with such small lattices one can get 2-state 
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signals even for crossovers, we use fourth-order 
Binder cumulants (-B4) 7 for tptp to get quanti- 
tative information about the transition: 
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where X is any observable. Since we use 5 noise 
vectors for each configuration to estimate ipip 
(and the number densities), we are able to ob- 
tain an unbiased estimate of B^ipip), as well as 
the susceptibility Xm- 



X X =V(X 2 -(X) 2 ). 



(3) 



The /3 of the transition is determined by the po- 
sitions of the maxima of the various susceptibili- 
ties, or from the minimum of B^ipip). (We have 
checked the consistency of these two methods.) 
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Figure 1. (3 C as functions of [x\ from minima of 
B4. The bottom point is for to = 0.03. The 
middle set of points and line are for m = 0.035. 
The top set of points and line are for to = 0.04. 



Figure shows the transition /3s, /3 C , as func- 
tions of fii for the 3 quark masses. The fits shown 
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which reflect the expected analyticity in fij. (3 C 
was obtained in each case from the position of 
the minimum of B^iptp). Ferrenberg-Swendsen 
reweighting was used to interpolate between the 
/3s of our simulations. These positions were in 
good agreement with those obtained from the 
maxima of the corresponding susceptibilities. 
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Figure 2. Binder cumulants for m = 0.035 as 
functions of /ij. The dashed line is at B4 = 1.604, 
the 3-d Ising value. 



For a lattice of infinite volume, the Binder cu- 
mulant B^ipip) is 3 for a crossover and 1 for a 
first-order transition. The critical endpoint is ex- 
pected to be in the universality class of the 3-d 
Ising model, for which B^-i/ji/j) = 1.604(1). While 
£?4 for a crossover or a first-order transition ap- 
proaches its infinite volume limit relatively slowly 
with volume, the approach for a second-order 
phase transition is fast and approaches the limit- 
ing value on relatively small lattices 8 . Figure [21 
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shows the Binder cumulant at the transition for 
our runs at m = 0.035 with hi — and hi = 0.2 
for both 8 3 x 4 and 12 3 x 4 lattices. Despite the 
sizable error bars, these are the results of runs of 
160,000 trajectories at the (3 closest to the tran- 
sition. 

-B4 at hi = lies above the Ising value while 
that at hi = 0.2 lies below this value for both lat- 
tice sizes. Straight lines through these 2 points 
for the two lattice sizes cross at some hi with 
0.0 < /ii < 0.2, and at a B4 value close to the 
Ising value. These preliminary results would indi- 
cate that there is a critical endpoint at hi = Vie 
with 0.0165 < $ E < 0.0231. 

For the m, fii values at which we have run and 
for j3s in the neighbourhood of the transition on 
an 8 3 x 4 lattice we estimate that the phase of the 
fermion determinant 9 obeys, 



9 2 ) « (5-12) M ?. 
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We have argued in that (cos 9) > 0.5 is a rea- 
sonable estimate of the range over which one can 
neglect the phase of the fermion determinant in 
determining the position and nature of the finite 
temperature transition l 5 . If so, the finite hi 
and 11 behaviour of this transition should agree 
for hi < 0.25 — 0.45, which includes our estimate 
for the critical endpoint for rn = 0.035. 

Using m c (0) = 0.0331 we then predict that 
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3. Conclusions 



We are simulating 3-flavour lattice QCD at fi- 
nite Hi and temperatures close to the finite tem- 
perature transition. 

We observe the slow decrease of (3 and hence 
temperature, of the finite temperature transition 
with fii. Arguing that, in this range of /ij, the 
dependence on hi an d 2/i should be the same we 
find a falloff about 20% slower than that observed 
by de Forcrand and Philipsen [3] At least some of 
the difference is due to omitting O(nj) terms. 

We have obtained preliminary evidence for a 
critical endpoint for m = 0.035 (just above 
m c (0)) from simulations at hi < 0.2. However, 



simulations on an 8 3 x 4 lattice at hi = 0.375 yield 
.84 = 1.64(2), which does not support this inter- 
pretation. This is probably because 0.375 is too 



close to He 



but this needs to be clarified. 



We have seen no evidence for a critical endpoint 
at m = 0.04. Whether this is consistent with 
what we see at m = 0.03 should also be checked. 

Our prediction for the hi dependence olm c (Hi) 
appears inconsistent with both de Forcrand and 
Philipsen [3| and the Bielefeld-Swansea collabo- 
ration [2], lying in between their inconsistent re- 
sults. 

Clearly much more needs to be done. 
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